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RESUMEN 

Se presenta la Teona de Poisson orientada al calculo de la Peligrosidad Sismica 
empleando Metodos Probabilisticos, con el proposito de que se conozca mas sobre este 
modelo y a futuro se utilicen con mas propiedad programas para evaluar la amenaza sismica. 

Luego se indica en forma resumida las zonas fuentes encontradas para el Ecuador 
para sismos corticales (Aguiar, 2009) y el calculo de las ventanas de tiempo en las cuales se 
considera que la informacion sismica es estable, aplicando el Metodo de Stepp en forma grafica 
y analltica. 

Posteriormente, se presenta las relaciones de recurrencia, para tres zonas fuente tipo 
de Ecuador que estan asociadas a sismos: corticales, interplaca tipo thrust e intraplaca. Para 
cada una de estas zonas fuentes se determina las relaciones de recurrencia de: Richter y 
Gutenberg (1944), Cornell y Vanmarcke (1969), Wesnousky (1994), Zoller et al. (2008). 

Finalmente se comparan las relaciones de recurrencia encontradas en este estudio 
para las tres zonas fuentes anotadas en el parrafo anterior, con las encontradas para Chile por 
Silva (2008) encontrando que las pendientes de las curvas de recurrencia son muy parecidas. 

Palabras Claves.- Relacion de recurrencia para Ecuador; Sismos interplaca tipo thrust; Sismos 
corticales; Sismos intraplaca. 


ABSTRACT 

Poisson Theory is presented focused on the calculation of Seismic Hazard using 
Probabilistic Methods, in order to know more about this model and to use programs with correct 
criteria to evaluate the seismic threaten. 

After that the source zones, found for Ecuador, are briefly indicated for cortical 
earthquakes and the calculation of time windows, where the seismic information is considered 
stable, applying Stepp Method graphically and analytically. 

Subsequently, the recurrence relations are presented, for three typical source zones of 
Ecuador which are related to: cortical earthquakes; interplate earthquakes of thrust type and 
interplate earthquakes. For each source, the recurrence relations are determined from: Richter 
and Gutenberg (1944), Cornell and Vanmarcke (1969), Wesnousky (1994), Zoller et al. (2008). 

Finally the recurrence relations, found in this study for the three source zones written 
down above, are compared with those found for Chile by Silva (2008), discovering that slopes of 
the recurrence curves are quite similar. 


Key Words.- Recurrence Relation for Ecuador; Interplate Earthquakes of Thrust type; Cortical 
Earthquakes; Interplate Earthquakes. 
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1. INTRODUCTION 

La sismicidad en una determinada area fuente puede ser modelada como un Proceso 
Estocastico de Poisson. Por este motivo se recuerda que el Proceso de Poisson, tambien 
conocido como Ley de los sucesos raros, es un proceso estocastico de tiempo continuo que 
consiste en contar eventos raros que ocurren a lo largo del tiempo; en este caso los eventos 
raros son los sismos que pueden tener diferentes magnitudes y pueden ocurrir en cualquier 
instante de tiempo. 

Se dice que el modelo de Poisson es un modelo sin memoria de tal forma que la 
ocurrencia de un sismo el 1 de enero de 2011, en una determinada area fuente no tiene ningun 
efecto sobre el numero de sismos que pueden ocurrir durante todo el mes de enero de 201 1 . El 
numero total de sismos que se registren en el mes es aleatorio y este valor mensual se 
considera constante durante todo el ano. En los modelos con memoria si influye el hecho de 
que se registro un sismo el 1 de enero de 2001, en el calculo del numero total de sismos 
mensual y anual. Por consiguiente en el modelo de Poisson, no se altera la evaluacion de 
la Peligrosidad Sismica, si se tiene un largo periodo de tiempo sin sismos. 

La funcion de distribucion del modelo de Poisson, viene definido de la siguiente 

manera: 


P(n > M)= — — - — (1) 

n\ 

Donde P(n > m) es la probabilidad de que en una determinada area fuente ocurran 

n sismos con magnitud mayor o igual a M en el intervalo de tiempo (0,?); A(m) es el 

numero medio de temblores con magnitud mayor o igual a M , que pueden registrarse en esa 
fuente. Un proceso de Poisson debe satisfacer las siguientes hipotesis. (Ruiz, 1985). 

i. Estacionariedad.- La probabilidad de un evento en un intervalo corto de tiempo 

{t,t + At) es aproximadamente A At para cualquier t. Es decir no importa que 

tiempo t se elige para hacer la estimacion. 

ii. No multiplicidad.- La probabilidad de 2 o mas eventos en un intervalo corto de tiempo 
es despreciable comparando con A At . 

iii. Independence.- El numero de eventos en cualquier intervalo de tiempo es 
independiente de el numero en cualquier otro intervalo de tiempo. 

La distribucion de Poisson es una funcion asimetrica en la cual el valor medio es igual a 

A y la desviacion estandar cr = Va .Si n = 0 (Probabilidad de no ocurrencia de sismos) se 
tiene, al reemplazar en la ecuacion (1), lo siguiente: 

P(0) = e~ At (2) 

Donde P( 0) es la probabilidad de no tener sismos de magnitud mayor o igual a M , 
en el intervalo de tiempo t . Por lo tanto, la probabilidad de que ocurran sismos con magnitud 
mayor o igual a M en el intervalo de tiempo t , vale: 


-At 


P{n) =\—e 


(3) 
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Tanto en la ecuacion (2) como en la ecuacion (3) se debio colocar X M en lugar de A , 
debido a que se esta encontrando la probabilidad de no tener o de tener sismos con magnitud 
mayor o igual a M . 

Sea T la variable aleatoria del tiempo en que ocurre un sismo de magnitud M . 
Entonces la probabilidad de que T exceda algun valor de t es igual a la que no ocurran sismos 
en el intervalo de tiempo t,oseaes igual a l-P T (t) = P(o). De donde: 

P T (t) = \-vop{-A, M t) t> 0 (4) 

Al derivar P T (t) con respecto al tiempo se halla la funcion de densidad de probabilidad 
f T (t ) , que vale: 


// (0 = dPl , (t) = K ex P (“ ( 5 ) 

at 

Donde t es el tiempo entre sismos con magnitud mayor o igual que M y /l M es la tasa 
de excedencia de la magnitud M . Luego el periodo de retorno T R (tiempo esperado entre 
sismos de magnitud M ) vale: 


T b = (6) 

A M 

2. RELACION DE RECURRENCIA 

La relacion de recurrencia de la actividad sismica, propuesta independientemente por 
Ishimoto-lda en 1939 y Richter-Gutenberg en 1944. Tiene la siguiente forma. (Richter, 1958) 

log N(M) = a-bM (7) 

Donde N(M)es el numero de sismos con magnitud igual o mayor que M ; a,b son 
parametros que se obtienen en cada area fuente aplicando regresion lineal. Se puede decir que 
a es igual al logaritmo del numero de sismos de magnitud mayor o igual que cero y b es la 
proporcion de sismos de una cierta magnitud. (Aguilar, 2001). De la ecuacion (7) se halla: 

N(M) = H) ia ~ bM} 
ln(N(M )) = (a-bM)* ln(l0) 

In (n(m )) = a ln(10) - b ln(10) M 

Se denomina: 

a = a ln(10) 
p = b ln(10) 

Con lo que se halla: 

\n(N(M)) = a-p M 

Sea A 0 el numero de sismos de magnitud mayor o igual a M 0 generados en promedio 
en una determinada area fuente en un tiempo t y sea A(M) el numero el numero de sismos de 


( 8 ) 

(9) 
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magnitud mayor a M en la misma fuente y en el mismo tiempo. Ahora al utilizar la ecuacion (9) 
se tiene: 

ln(/l 0 ) = a - ft M 0 
ln(/l(M)) = a - P M 

Al restar estas ecuaciones se encuentra: 


ln(/l 0 ) — ln(/l(M)) = ~J3M 0 +J3M 


Al cambiar de signo y facturando el lado derecho se tiene: 


De donde: 


In 


\ K j 




A(M) = A 0 e~ /3(M ~ M ° ) (10) 

El problema que tiene la ecuacion (10) es que no tiene no tiene ningun Ifmite superior 
para la magnitud maxima M v . 

Cornell y Vanmarcke (1969) presentan un modelo mas realista para encontrar A(M) en 
lugar de la ecuacion (1 0) y es el siguiente: 


A(M) = A 0 


e -PM _ e -PM v 
- e ~ pMv 


( 11 ) 


La ecuacion (1 1) se utiliza en algunos programas que existen para evaluar Peligrosidad 
sismica, como es el caso del Programa CRISIS, en sus distintas versiones. Ordaz et al. (2007), 
cuando se utiliza la funcion de distribution de Poisson. Si se dibuja la ecuacion (11) en papel 
semilogaritmico se vera que para magnitudes pequenas es una recta con una pendiente igual a 
— P lo que esta en concordancia con lo observado por Gutenberg y Richter (1944), para 

magnitudes altas la recta se transforma en una curva concava hacia abajo y A(m) = 0 para 
magnitudes mayores a M v . 


El calculo de los parametros A 0 , /3, M v se lo realiza en base a los datos del catalogo 

sfsmico pero teniendo en cuenta el problema de la completez de la information sismica que se 
presenta en el siguiente apartado. Pero al margen de esto es importante tener en cuenta las 
siguientes observaciones. Ordaz (2002). 

• El tiempo de observation de sismos instrumentales es muy limitado. En el caso de 
Ecuador, las redes sismicas empezaron a funcionar en el periodo 1980-2000. Ademas 
que existe incertidumbre en la ubicacion de los epicentros, profundidad focal y magnitud 
con respecto a lo que proporciona el USGS, para algunos eventos sfsmicos. 

• Los catalogos sfsmicos vienen en varias escalas de magnitud; M L para los sismos 

registrados antes de 1960; M b , M s para los sismos registrados entre 1960 y el 2000; 
M w para los sismos posteriores al 2000. 

Por todo esto es conveniente comparar los parametros indicados, con los encontrados 
en otras regiones con tectonica similar y tomar una decision basada en la experiencia de quien 
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realiza el estudio de peligrosidad sismica. Algo similar sucede cuando se calcula un edificio, los 
resultados que reportan los programas son referenciales, quien toma la decision de como armar 
la estructura es el Proyectista Estructural basado en su experiencia y conocimientos, si no fuera 
de esta manera, estarla demas el Proyectista ya que existen programas que dibujan los pianos 
estructurales pero esto no es asl. 



3. COMPLETEZ DE LA INFORMACION SISMICA 

Si se observan los catalogos slsmicos se vera que muy pocos sismos estan registrados 
en el perlodo 1900-1960 de magnitudes comprendidas entre 4 y 4.5, por decir algo, una 
cantidad mayor se tienen registrados entre 1960-1990 y bastantes son los sismos registrados 
en el rango de las magnitudes anotadas con posterioridad a 1990. A esta falta de informacion 
en ciertas ventanas de tiempo se conoce con el nombre de completez de la informacion 
sismica. 

Stepp (1973) propuso un procedimiento para determinar ventanas de tiempo en los 
cuales la informacion sismica puede considerarse completa, basado en el criterio de estabilidad 
de la tasa de recurrencia, para ciertos rangos de magnitud. 
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El procedimiento propuesto por Stepp ha sido utilizado en trabajos recientes para 
evaluar la peligrosidad sismica en: Chile, (Silva, 2008) Centro America (Molina et al. 2008). En 
el ultimo trabajo se aplico el Metodo de Stepp en forma grafica. En Aguiar et al. (2011) estan 
detallados la forma de calculo de las ventanas de tiempo en forma grafica y en forma analftica, 
aplicados a las zonas fuentes del Ecuador. 

Con el proposito de explicar estos dos metodos se presenta el calculo de las ventanas 
de tiempo para la Zona Fuente 2, indicada en la figura 1, destacando que no es objetivo de este 
artfculo describir como se obtuvieron las zonas fuentes ya que estan indicadas en Aguiar y 
Castro (2009) una primera version reciente y en Aguiar (2010) una segunda version. 


COMPLETEZ ZONA 2 



Ano 

Figura 2 Numero acumulado de sismos normalizados en la fuente cortical 1 . 

En la figura 2 se presenta el numero acumulado de sismos, por rangos de magnitud, 
desde 1 540 hasta el 2000 pero normalizados a la unidad. La ventana de tiempo es aquella en la 
cual la pendiente es constante, esto es lo que se debe observar en cada rango de magnitud; en 
la figura 2 se aprecian varias pendientes con lo cual queda de manifiesto el problema de la falta 
de datos en ciertas ventanas de tiempo. 

En la figura 3 se ha dibujado los sismos acumulados y normalizados, con magnitudes 
entre 4.5 y 5.0, todo esto es del area fuente 2. Se puede ver en esta grafica que a partir de 1990 
se tiene una sola pendiente Lugo la ventana de tiempo para los sismos entre las magnitudes 
indicadas es de 20 anos. Lo propio se hace para cada rango de magnitud. 

Ahora, en el Metodo de Stepp, propiamente dicho se encuentra la desviacion estandar 
para diferentes ventanas de tiempo y para cada rango de magnitud. En la parte superior de la 
figura 4 se presenta este calculo para sismos con magnitudes entre 4.5 y 5.0 de la fuente 2. 

Despues se debe dibujar una recta paralela a — que pase por esos puntos. Donde 

T es el tiempo. El punto en el cual se logra el mejor ajuste es aquel que determina el ancho de 
la ventana de tiempo que en este caso es 1980. Se esta trabajando en escala log-log. 
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COMPLETEZ ZONA 2 Magnitud (4.5 - 5.0) 



Figura 3 Resultado del estudio de completez para sismos con magnitudes entre 4.5 y 5.0 


COMPLETEZ STEPP ZON A 2 Magnitud (4.5 - 5.0) 
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Figura 4 Resultados con el Metodo de Stepp, para sismos con magnitud 4.5 y 5.0 

El Metodo de Stepp considera completa la muestra en cada rango de magnitud y 
especfficamente en una determinada ventana de tiempo, cuando el valor medio R t es 

1 

constante y no depende del tiempo en consecuencia se comporta de la forma —j = . El valor de 

& 

R, se obtiene dividiendo el numero total de sismos encontrados en una determinada ventana 
de tiempo para el ancho de la ventana de tiempo. 

En la tabla 1 se indican los resultados hallados para la zona fuente 2, del estudio de 
completez de la informacion sismica. Aguiar etal. (2011). 
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Tabla 1 Ventanas de tiempo halladas para los sismos de la Zona Fuente 2. 


MAGNITUD 

4. 0-4. 5 

4. 5-5.0 

5.0-5. 5 

5.5-6. 0 

6. 0-6. 5 

6.5-7.0 

7. 0-7. 5 

7. 5-8.0 

VENTANA 

DE 

TIEMPO 

(anos) 

20 

20 

20 

30 






Las ventanas de tiempo que se han obtenido son muy cortas, con lo que se incrementa 
la incertidumbre en la evaluacion de la peligrosidad slsmica ya que por ejemplo, para el caso de 
la tabla 1 se considera en el mejor de los casos 30 anos de datos sfsmicos y se deja de lado 
valiosa informacion de sismos registrados en el pasado de gran magnitud. 

4. RELACIONES DE RECURRENCIA 

Con el proposito de ilustrar la forma de calculo de las relaciones de recurrencia; en la 
tabla 2 se presenta el calculo de los parametros a, b , de la ecuacion (7), con los datos 

slsmicos de la zona fuente 2 hasta diciembre de 2010, para el ajuste se considera que la 
Magnitud es la variable X y que el numero de sismos de un determinado rango de magnitud 
dividido para la longitud de la ventana de tiempo es la variable Y. 


Tabla 2 Calculo de los parametros a, b de la ley de recurrencia de Richter. 


N 

Magnitud 

M w 

Ventanas 
T (anos) 

Numero 
de sismos 
(NS) 

Y=NS/T 

Log Y 

X*log Y 

X 2 

1 

4. 0-4.5 

4.0 

20 

158 

7.900 

0.898 

3.591 

16.00 

2 

4. 5-5.0 

4.5 

20 

25 

1.250 

0.097 

0.436 

20.25 

3 

5. 0-5.5 

5.0 

20 

6 

0.300 

-0.523 

-2.614 

25.00 

4 

5. 5-6.0 

5.5 

30 

2 

0.067 

-1.176 

-6.469 

30.25 

Sumatoria 

19 


191 

9.517 

-0.533 

-5.056 

91.50 


Las ecuaciones con las cuales se obtienen los parametros a, b , son las siguientes 
(Ajuste de rectas por mmimos cuadrados). 


^log y = N g + 
^xlogy = + £> 


( 12 ) 


Al reemplazar los datos de la tabla 2, en la ecuacion (12) se halla: a = 6.3228 
& = 1.3682 Ahora bien si se grafica la ecuacion (7) con los valores a, b encontrados se 
obtendra una recta. 

La magnitud maxima registrada en una fuente, no necesariamente se encuentra dentro 
de la ventana de tiempo considerada o puede ser que la magnitud maxima esperada en la 
fuente sea mayor a la magnitud maxima registrada. En este caso hay dos planteamientos el 
primero que considera que la ley de Gutenberg-Richter (GR) es bien apropiado para sismos de 
magnitud moderada y baja pero para sismos de magnitud alta se plantea un escalamiento 
prolongando la recta de la relacion de GR hasta la intersection con la maxima magnitud 
esperada. Zoller etal. (2008). 

El segundo criterio a partir del final de la recta que se obtiene con la ecuacion (7) se 
traza una linea horizontal hasta el punto de code, con una recta perpendicular a la magnitud 
maxima. Wesnousky (1994). Que en definitiva es otra forma de reconocer que GR es adecuado 


Relaciones de recurrencia, aplicacion a la peligrosidad si'smica del Ecuador 


209 


para magnitudes medias y bajas y la forma de escalamiento para magnitudes altas es con la 
recta paralela, ver figura 5. 

Para aplicar el criterio de Cornell y Vanmarcke (1969) indicado en la ecuacion (11) se 
debe encontrar el valor de A 0 que es el numero de sismos con magnitud mayor a 4.0, al ano. 
Para ello se debe obtener: 


a = a *ln(10) = 14.5588; /? = £>*ln(10) = 3.1504; A 0 = exp (a — ft M 0 ) = 7.0800 


En la figura 5 se presentan las relaciones de recurrencia aplicando los tres criterios, a la 
izquierda el de Cornell y Vanmarcke (1969), la magnitud maxima registrada en la fuente 2 es 
6.26.; al centra se tiene la relation de recurrencia con el criterio de Wesnousky (1994) y a la 
derecha con el criterio de Zoller et al. (2008). De los tres el criterio de Cornell y Vanmarcke 
(1969) es el que se utiliza. 



Wesnousky (1994) ZONA 2 



Zoller et al (2008) ZONA 2 


a) Cornell y Vanmarcke (1969) 


b) Wesnousky (1994) 



Figura 5 Relaciones de recurrencia, para zona fuente 2. 


Si se reemplaza en la ecuacion (1) A 0 = 7.08 , t = 1 y se encuentra la probabilidad de 

ocurrencia anual de sismos con una distribucion Poisson, se obtiene la curva indicada en la 
figura 6, que representa la probabilidad de ocurrencia de sismos de magnitud mayor o igual a 4 
en la fuente 2. 



Figura 6 Probabilidad anual de ocurrencia de sismos de magnitud mayor a 4 en fuente 2. 


Para A 0 «7se obtiene la mayor probabilidad que es de aproximadamente 0.15; este 

es el valor medio de la funcion de distribucion que no es simetrica. La probabilidad de que en la 
fuente 2 no ocurra ningun sismo de magnitud mayor a 4 es practicamente cero (esto se observa 

del grafico); ahora, al utilizar la ecuacion (2) se halla que la probabilidad vale 8.41*10^. 
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Tambien se aprecia que la probabilidad de tener mas de 16 sismos al ano de magnitud mayor a 
4 es nula. 



Figura 7 Relaciones de recurrencia para las zonas fuentes 6, 14 y 16. 

Para posteriormente comparar las relaciones de recurrencia de los sismos de Ecuador 
con los sismos de Chile, se indican los valores hallados para Chile. Silva (2008). 

log (iv) = 5.09 — 0.865 M s 

log (iv) = 5.39 — 0.965 M s (13) 

log (N) = 5.08-1.061 M s 

La primera de las ecuaciones identificadas con el numero (13) es para sismos interplaca 
tipo thrust; la segunda para sismos intraplaca y la tercera para fallamiento cortical. Se destaca 
que estan trabajando con magnitud M s y tambien que para magnitudes altas se tiene que 

M w ~ M s . En el Ecuador se trabajo con magnitudes M w . 
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5. RESULTADOS 

La figura 7 presenta relaciones de recurrencia similares a los presentados en la figura 5 
pero para las zonas fuentes 6, 14 y 16. Los comentarios que se pueden hacer a estas graficas 
son similares a los ya realizados, cuando se presento los resultados de la fuente 2. Unicamente 
quedaria pendiente indicar que la zona fuente 14 corresponde a sismos interplaca tipo thrust y 
la zona fuente 16 a sismos intraplaca asociados con la subduccion. En la figura 1 estan las 
fuentes asociadas a sismos corticales. 

Ahora bien, en la tabla 3 se presentan los valores a, b de la ley de recurrencia de 
Richter; el valor de A 0 que es la tasa de ocurrencia anual de sismos con magnitud mayor o 

igual que 4; el parametro J3 = b ln(10) y la magnitud maxima M MAX que ha sido registrada en 
la zona fuente. 


Tabla 3 Parametros que definen la sismicidad en las zonas fuentes 2, 6, 14 y 16. 


ZONA 

a 

b 

A 0 

P 

M max 

2 

6.3228 

1.3682 

7.0800 

3.1504 

6.26 

6 

7.8123 

1 .5444 

43.1030 

3.5562 

6.38 

14 

4.8273 

0.9304 

12.7558 

2.1423 

7.10 

16 

4.8112 

0.9499 

10.2669 

2.1873 

7.50 


Las fuentes 2 y 6 estan asociadas con fallamiento cortical, la 14 con interplaca tipo thrust 
y la 16 con intraplaca de subduccion. Con este recordatorio en la figura 8 se presentan las 
relaciones de recurrencia para las fuentes indicadas, aplicando el criterio de Zoller et al. (2008) 
mas con la idea de tener Ifneas rectas. En esta figura se observa que las fuentes 14 y 16 tienen 
rectas que son practicamente paralelas (valor similar de b)\ la fuente 14 tiene mayor 
sismicidad que la 16 (valor de a es mayor en la fuente 14 que en la 16). 



En base a la figura 8 se puede indicar que en el Ecuador la mayor cantidad de sismos 
de magnitudes medias y altas corresponden a sismos interplaca tipo thrust; seguidos por los 
sismos intraplaca de subduccion y finalmente los sismos corticales. 
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Figura 9 Relaciones de recurrencia para Ecuador y Chile. 


En la figura 9 se comparan las relaciones de recurrencia encontradas para Ecuador en 
las zonas fuente 2, 14 y 16, con las encontradas para Chile para sismos corticales, interplaca 
tipo thrust e intraplaca. En forma general se puede decir que las dos rectas que contiene cada 
grafica son paralelas y que Chile tiene una mayor actividad sfsmica que el Ecuador. Pero lo 
interesante es ver la forma de las rectas que son parecidas lo cual significa que los valores 
encontrados para el Ecuador son correctos. 

6. CONCLUSIONES 

Se ha presentado el marco teorico con un ejemplo, de la funcion de distribution de 
Poisson orientada al calculo de la Peligrosidad sfsmica empleando Metodos Probabilfsticos. 
Luego se ha indicado en forma resumida las zonas fuentes encontradas para el Ecuador para 
sismos corticales y el calculo de las ventanas de tiempo en las cuales se considera que la 
informacion sfsmica es estable, en cada rango de magnitud aplicando el Metodo de Stepp en 
forma grafica y analftica. Todo esto con el proposito de que el lector pueda seguir el calculo de 
las relaciones de recurrencia aplicando las propuestas de: Richter y Gutenberg (1944); Cornell y 
Vanmarcke (1969); Wesnousky (1994); Zoller et al. (2008), que es el objetivo del artfculo. 

Posteriormente, se presentan las relaciones de recurrencia, para tres zonas fuente tipo 
de Ecuador que estan asociadas a sismos: corticales (zona fuente 2); interplaca tipo thrust 
(zona 14) e intraplaca (zona 16). Para cada una de estas zonas fuentes se determina las 
relaciones de recurrencia con cada uno de los modelos indicados en el parrafo anterior. 
Unicamente para el modelo de Richter y Gutenberg se han comparado los resultados obtenidos 
con los hallados por Silva (2008) en Chile. Del estudio realizado se desprenden las siguientes 
conclusiones: 
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• Es muy importante que las relaciones de recurrencia tengan una magnitud maxima 
esperada. En este contexto los modelos de: Cornell y Vanmarcke (1969); Wesnousky 
(1994); Zoller etal. (2008), son adecuados. 

• En el Ecuador las relaciones de recurrencia para sismos interplaca tipo thrust e 
intraplaca presentan pendientes paralelas. 

• Se nota que las mayor sismicidad se registra en sismos interplaca tipo thrust, seguido 
de los intraplaca y finalmente los sismos corticales. Esta tendencia es similar a la 
hallada en Chile por Silva (2008). 

• Las pendientes de las relaciones de recurrencia halladas en el Ecuador para sismos 
interplaca de bajo angulo de buzamiento e intraplaca tienen pendientes muy parecidas 
a las encontradas en Chile. 
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